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Abstract 

We propose a 4-dimensional Kaluza-Klein approach to general relativity in 

^^ , the (2,2)-splitting of space-time using the double null gauge. The associated 

Q"^ i Lagrangian density, implemented with the auxiliary equations associated with 

qh| the double null gauge, is equivalent to the Einstein-Hilbert Lagrangian density, 

}^, since it yields the same field equations as the E-H Lagrangian density does. It 
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is describable as a (l-l-l)-dimensional Yang-Mills type gauge theory coupled to 
(l-l-l)-dimensional matter fields, where the minimal coupling associated with 
the infinite dimensional diffeomorphism group of the 2-dimensional spacelike 
fibre space automatically appears. The physical degrees of freedom of gravi- 
tational field show up as a (l-l-l)-dimensional non-linear sigma model in our 
Lagrangian density. Written in the first-order formalism, our Lagrangian den- 
sity directly yields a non-zero local Hamiltonian density, where the associated 
time function is the retarded time. From this Hamiltonian density, we obtain 
a positive-definite local gravitational energy density. In the asymptotically 
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flat space-times, the volume integrals of the proposed local gravitational en- 
ergy density over suitable 3-dimensional hypersurfaces correctly reproduce 
the Bondi mass and the ADM mass expressed as surface integrals at null 
and spatial infinity, respectively, supporting our proposal. We also obtain the 
Bondi mass-loss formula as a negative-definite flux integral of a bilinear in the 

gravitational currents at null infinity. 
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I. INTRODUCTION 

The exact correspondence of the Euchdean self-dual Einstein's equations to the equa- 
tions of motion of 2- dimensional non-linear sigma models with the target space as the area- 
preserving diffeomorphism of 2-surface 0,^ has inspired us to look further into the intriguing 
question whether the full-fledged general relativity of the 4-dimensional space-time can be 
also formulated as a certain (l+l)-dimensional field theory. Recently we have shown that 
such a description is indeed possible, and constructed the action principle [§|-§1 in the frame- 
work of the 4-dimensional Kaluza-Klein theory in the (2,2)-splitting. In this approach, the 
4-dimensional space-time, at least for a finite range of space-time, is viewed as a fibred man- 
ifold that consists of the (l+l)-dimensional "space-time" and the 2-dimensional "auxiliary" 
fibre space. 

There are certain advantages of this 4-dimensional KK approach to general relativity 
in the above splitting, which led us to develop this formalism. We list a few of them. 
First of all, in (H-l)-dimensions, there exist a number of field theoretic methods recently 
developed thanks to the string-related theories. Hopefully the rich mathematical methods in 
(l+l)-dimensions might also prove useful in studying general relativity in the (2,2)-splitting, 
classical and quantum. Moreover, in this KK formulation, general relativity can be viewed as 
a (l+l)-dimensional gauge theory with the prescribed interactions and auxiliary equations. 
Since the major advantage of gauge theory formulation is that gauge invariant quantities 
automatically solve Gauss-law equations associated with the gauge invariance, the problem 
of solving constraints of general relativity could be made even trivial, at least for some of 
them. Furthermore, this formulation allows us to forget about the space-time picture of 
general relativity; instead, it enables us to study general relativity much the way as we do 
for Yang-Mills theories coupled to matter fields in (H-l)-dimensions, putting the space-time 
physics into a new perspective. 

This (l+l)-dimensional method, however, is not entirely new since it was virtually used 
in analyzing gravitational waves W\M, and was further developed in the spin-coefficient 



formalism PH12[ and the null hypersurface formalism |r^-[T9| . In these formalisms, a special 



gauge which we may call the double null gauge is chosen such that two real dual null 
vector fields whose congruences span the (l+l)-dimensional submanifold are singled out, 
and the Einstein's equations are spelled out in that gauge. A characteristic feature of these 
formalisms, among others, lies in that the true physical degrees of freedom of gravitational 
field show up in the conformal 2-geometry of the transverse 2-surface [pip!3 HT7|j20|l . This 



feature, that has been particularly useful for studying the propagation of gravitational waves 
in the asymptotically flat space-times, further motivated the canonical analysis [p!8iplH 26[| 
of the null hypersurface formalism, in the hope of getting quantum theory of gravity by 
quantizing the true physical degrees of freedom of gravitational field. 

In view of these advantages of the KK formalism and the null hypersurface formalism, it 
therefore seems worth combining both formalisms to see what could be learnt more about 
general relativity. In this article, we shall present such a formalism. In this approach, 
it is the (l+l)-dimensional submanifold spanned by two real dual null vector fields that 
we imagine as "space-time" and the remaining transverse 2-surface as the "auxiliary" fibre 
space. As a by-product of our KK approach in the double null gauge, we obtain a new result 
which we report in this article. Namely, we propose a positive- definite local gravitational 
energy density in general relativity without referring to the boundary conditions, and show 
that the volume integrals of the proposed energy density over suitably chosen 3-dimensional 
hypersurfaces correctly yield the Bondi and the ADM surface integral at null and spatial 
infinity in the asymptotically fiat space-times, respectively. We also obtain the Bondi mass- 
loss formula as a negative-definite fiux integral of the gravitational degrees of freedom at 



null infinity |7|,|8|JTT[| . The proposed local gravitational energy density comes directly from 
the local Hamiltonian density of general relativity described as the 4-dimensional KK theory 
in the double null gauge. The associated time function is the retarded time, and has the 



physical interpretation |^ as the phase of the local gravitational radiation in situations 
where gravitational waves are present. 

This article is organized as follows. In section II, we present the 4-dimensional KK theory 



in the (2,2)-splitting, using the double null gauge. It will be seen that, even when the gauge 
symmetry is an infinite dimensional symmetry such as the group of diffeomorphisms, the 
KK idea is still useful by showing that the KK variables transform properly as gauge fields 
and tensor fields under the corresponding gauge transformations p|-|6|,p8|. Next, we present 
the Lagrangian density for general relativity in the double null gauge, which is implemented 
by the auxiliary equations associated with the double null gauge using the Lagrange mul- 
tipliers. This Lagrangian density is equivalent to the Einstein-Hilbert Lagrangian density, 
since it yields the same field equations as the E-H Lagrangian density does. We shall present 
the 10 Einstein's equations in the double null gauge. Moreover, we shall find that 2 of these 
10 equations are in a form of the Schrodinger equation, reminiscent of the Brill wave equa- 
tion [p9|-|31|, i.e. the initial value equation of the axi-symmetric gravitational waves at the 
moment of time symmetry. 

In section III, we shall present this Lagrangian density in the first-order formalism, 
with the retarded time identified as our clock variable. This immediately leads to the local 
Hamiltonian density and thus to the local gravitational energy density that we are interested 
in. The proposed local gravitational energy density is positive-definite. We shall further show 
that the volume integral of the proposed local gravitational energy density over a suitably 
chosen 3-dimensional hypersurface become a surface integral, using the vacuum Einstein's 
equations and the Bianchi identities. In the asymptotically fiat space-times, this surface 
integral becomes the Bondi and the ADM surface integral defined at null and spatial infinity, 
respectively. We also derive the Bondi mass-loss formula from the proposed gravitational 
energy density. 

In Appendix A, we shall describe the general (2,2)-splitting of space-time, and present 
the resulting E-H Lagrangian density without picking up a special gauge, as we need it when 
we wish to obtain the field equations from the variational principle. In Appendix B, we shall 
introduce the covariant null tetrads, as we shall use them in Appendix C where we show 
that the proposed volume integrals in section III can be expressed as surface integrals. In 
Appendix D, we shall show that /tj_, which appears when we discuss the Bondi mass- loss, is 



positive-definite. 

II. THE LAGRANGIAN DENSITY IN THE DOUBLE NULL GAUGE 

In this section we combine the null hypersurface formalism with the 4-dimensional 
Kaluza-Klein approach where space-time is viewed as a fibred manifold, i.e. a local product of 
the (l+l)-dimensional base manifold and the 2-dimensional fibre space. Let the vector fields 
djdX^ = {d/du, d/dv, d/dy"") (a = 2, 3) span the 4-dimensional space-time. In a Lorentzian 
space-time we consider here, there always exist two real null vector fields, which we may 
choose orthogonal to the 2-dimensional spacelike surface N2 spanned by d/dy"'. Following 
the KK idea ||3^ , the two null vector fields can be represented as the linear combinations of 



these basis vector fields 

B B B B 

^-A^i-, and ^-A_^^, (2.1) 
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for some functions A^{u,v,y). Since these null vector fields are assumed to be normal to 
N2, the line element may be written in a manifestly symmetric way as follows; 

ds^ = -2dudv + (t)ab{A^du + A^dv + dy''){A_^du + A^dv + dy''), (2.2) 

where 0ab(w, v, y) is the 2-dimensional metric on N2. Notice that, as a consequence of picking 
up two null vector fields normal to N2, 2 out of the 10 metric coefficient functions were 
gauged away in (|2.2| ). In addition, one more function was removed from (|2.2| ) by choosing 
the coordinate v such that Cdv' = dv for some function C, i.e. by choosing C = 1. The 
elimination of these 3 functions may be viewed as a partial gauge-fixing of the space-time 
diffeomorphism, and may be better understood in terms of the dual metric, which we may 
write 

g++=g- = 0, g+-=g-+ = -l, g+" = A_^, 

g-" = A_^, g"^ = (f' - 2A^At (2.3) 



That g^^ = g = means that du and dv are dual null vector fields, and that g~^~ = — 1 
is a normalization condition for v, given an arbitrary function u. We shall call this gauge 
as the double null gaugeQ, and general relativity formulated in this gauge is referred to as 
the double null formalism |r^ -p!7| . The 3-dimensional hypersurface defined hj u = constant 



is a null hypersurface since it is metrically degenerate; within each null hypersurface the 
2-dimensional spacelike space N2 defined hj v = constant is transverse to both du and dv. 
In order to see whether this 4-dimensional KK program is justifiable in the absence of any 
Killing symmetry, as is the case here, we have to first examine the transformation properties 
of (pab and A^ in ( |2.2|) under the action of some group of transformations associated with 



N2. The most natural group of transformations associated with N2 is the diffeomorphisms 
of N2, i.e. diffA^2- Under the diffiV2 transformation 

y'' = y'\u,v,y), u' = u, v' = v, (2.4) 

these fields must transform as 

dy'^ dy 
'P'abiu, V, y') = ^-^^0cd(M, V, y), (2.5a) 

Ai%u, V, y') = ^^± («, V, y) - d^y'\ (2.5b) 

so that the line element ds^ remains invariant [§-0]. Under the corresponding infinitesimal 
t r ansf ormation 

5y'^ = i'^{u,v,y), 5u = 5v = Q, (2.6) 

where 1^" is an arbitrary function, we find 

^0a6 = -[^,0]a6, (2.7a) 

8A^ = -D^C = -d±e + [A±, e]^ (2.7b) 



^We notice that this double null gauge is valid only for a finite range of space-time. See for 



instance |1S,14,30|. 



where A± := A_^da and ^ := ^""da- Here the brackets are the Lie derivatives associated with 
diffiVs, 

[e, <PU = i"dcct>ab + {daO(t>cb + {d,Oct>ac, (2.8a) 

[A^,ir = Ald,C-CdcA^. (2.8b) 



This observation tells us two things. First, diffA'"2 is the residual symmetry |Tl| of the line 
element (|2.2| ) which survives even after the double null gauge was chosen. Second, diffiV2 
should be viewed as a localQ gauge symmetry of the Yang-Mills type, since A^ and (pah 
transform as a gauge field and a tensor field under the diffiV2 transformations, respectively. 
This feature is rather surprising, since the KK variables were often thought to be useful for 
higher dimensional gravity theories where the degrees of freedom associated with the extra 
dimensions are suppressed in one way or another^. In our 4-dimensional KK approach to 
general relativity, we leave all the "internal" degrees of freedom intact so that all the fields in 
( ^.2[ ) depend on all of the coordinates (u, f , y""). In spite of these generalities, the Lagrangian 
density associated with the metric (|2.2|) can be still identified, a la Kaluza-Klein, as a gauge 
theory Lagrangian density defined on the (l+l)-dimensional "space-time", with diffiV2 as 
the associated local gauge symmetry, as we shall see shortly. 



The metric (|2.2|) in the double null gauge can be obtained from the general KK line 
element 

ds"^ = -f^^dx^'dx" + (pabiAyx'' + dy''){Aj'dx'' + dy^), (2.9) 



where /x, z/ = 0, 1 and a,b = 2, 3. From this we obtain (|2.2| ) by introducing the retarded and 



advanced coordinate {u,v) and the fields A^, 



u = —{x^ -x^), t. = — (x° + x^), (2.10a) 

A^^ = ^{A,^tA,'^), (2.10b) 



^Here local means local in the (l+l)-diniensional "space-time". 
^See however [28,33|. 



assuming the (H-l)-diniensional "space-time" metric 7^jy to be 



1+- = -1, 7++ = 7— = 0, 



(2.11) 



in the (m, f )-coordinates, where +(— ) represents u{v). In Appendix A, the general E-H 
Lagrangian density 0-§[ for the metric ( |2.9| ) and the prescription how to obtain it are 
presented. Using the "ansatz" ( |2.1lD , we can easily show that the general E-H Lagrangian 
density (|A32|) reduces to the following expression corresponding to the metric (|2.2| ). If we 
neglect the auxihary equations associated with the double null gauge (|2.11|) for the moment, 
it is given by [§] 



Co 



^ >abF^^F^'_ + ]-r' (!>''' [{D+ct>ac){D-Ct>M) " (Z^+^afa) P-0ed) } 



(2.12) 



where we ignored the surface terms. Here = det 0ab, and F^_^ is the diffiV2-valued field 
strength, and D±(j)ab is the diffA^2-covariant derivative defined as 



D±(j)ab = d±(j)ab - [A±, 0]afe, 

where [^4+, A_Y and [A^, (f)\ab are the Lie derivatives defined as 
[A+, AJ]" = A^d^A^ - A_'d,A_^, 

[A±, (PU = A^d,(f)ab + {daA^)(f),b + {dbA^)(j)ac, 



(2.13a) 
(2.13b) 



(2.14a) 
(2.14b) 



respectively (see Appendix A). Recall that, in the double null formalism of general relativity, 
the transverse 2-metric with a unit determinant, i.e. the conformal 2-geometry, is the 
two physical degrees of freedom of gravitational field ||8|J13|-|T7| . Since we are essentially 
reformulating the double null formalism from the KK point of view in this article, we would 
like to see first of all what the Lagrangian density looks like when written in terms of 
the conformal 2-geometry. Let us therefore decompose the 2-metric (pab into the conformal 
classes 



■>ab 



flpab, {fl > and detpab = l) 



(2.15) 



where pab is the conformal 2-geometryQ of the transverse 2-surface N2- If we define a by 
a := lnr2, the second term in ( p.l2| ) becomes 

= -Q-\D+Q){D^n) + hlp'''p^\D+p,,){D^pM) 

= -e'^(Z}+a)(D_a) + ^e>V'(^+Pac)(/^-PM), (2.16) 

where D±VL, D±a, and D±pab are the diffiV2-covariant derivatives 

D^Vl = d±n - [A±, n], (2.17a) 

D±a = d±a-[A±,a], (2.17b) 

D±Pab = d±Pab - [A±, p]ab, (2.17c) 

and [A±,fi], [A±,cr], and [A±, p]ab are given by 

[A±, n] = A^da^ + {daA^)n, (2.18a) 

[A±,a] = A^daa + d,A^, (2.18b) 

[A±, p]ab = A^d.pab + {daA^)p,b + {dbA^)pac - {dcA^)pab, (2.18c) 

respectively. Here 9aA "-terms are included in the Lie derivatives, since Q and pab are tensor 
densities of weight —1 and +1 under diffA^2; respectively. Thus (|2.12|) becomes 



Co = -e^'^PabF^^F^' - e'^iD+a)iD_a) + -eVV'(^+Pac)(/^-P6d). (2.19) 

In order to find the correct variational principle that yields all the 10 Einstein's equations, 
however, we must implement (|2.19| ) with the auxiliary equations associated with the double 



^This may be viewed as a finite analogue of the physical transverse traceless degrees of freedom of 
the spin-2 fields propagating in the flat space-time [Q. The double null gauge may be also viewed 
as an analogue of the Coulomb gauge in Maxwell's theory, where the physical degrees of freedom 
are the transverse traceless vector potentials [p^ ]. 
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null gauge. These equations can be found by first varying the general E-H Lagrangian 
density ( |A32D with respect to 7^~*',7 , and 7"^~, and then plugging the double null gauge 
( p. lip into the resulting equations. They are found to be 

C± := Dla + ]^{D±af + -^p'^' p''\D±pac){D±pM) = 0, (2.20a) 

Co := -U'' PabF+^F^"^ + D+D_a + D_D+a + 2{D+a){D_a) + R^ 

= 0, (2.20b) 

respectively, where i?2 := <P°''^Rac is the scalar curvature of A^2- Remarkably, the equations 
C± = in ( p.20a| ) may be written in a form of the Schrodinger equation. Using the identity 

e^/^[Dla + -{D^af) = 2 Z}|e'^/^ (2.21) 

the equations C± = become 

Z)^e'^/2 + 4e'^/2 _ 0, where 4 - ^pV''(^±Pac)(/^±Pbd) > 0. (2.22) 

o 

That K^, a bilinear in the currents of the gravitational degrees of freedom, is positive- 



definite can be shown easily (see Appendix D). The equations ( |2.22[ ) are the analogues 



of the Brill wave equation^ p5|-pT|], as they are of the Schrodinger equation type for a 
wave function corresponding to a state of zero energy in the potential — k±, coupled to 
the external gauge fields A_^\ Thus, viewed as a scattering problem, the scattering data 
e"^'^ in (|2.22|) is an auxiliary field that can be determined by the potential — k^ up to 



^Our Schrodinger equations look like one-dimensional wave equations coupled to gauge fields, but 
actually they are 3-diniensional partial differential equations like the Brill wave equation, due to 



the Lie derivatives along A± = A^da in (2.22). The wave function in (2.22) is related to the 



confornial factor of the 2-diniensional wavefront, the spatial projection of the null hypersurface 
u = constant, rather than that of spacelike hypersurface. But it should be also mentioned that, 
for the metric ( |2.2| ), the area measure of the 2-dimensional wavefront and the volume measure of 
3-dimensional null hypersurface u = constant are the same [ pq| . 
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some integral "constant" functions. The generic behaviors of solutions of the 2 Einstein's 
equations C± = are therefore expected either of the scattering type, or of the bound-state 
or resonance type, corresponding to the asymptotically flat space-times or spatially closed 
universes, respectively, on a par with the Brill wave equation. 
The correct variational principle is now given by 

C = h^''pabF_,^F^' - e"{D+a){D^a) + ^eW'(^+Pac)(l^-P6d) + E ^"C-^, (2.23) 

where A"'s are the Lagrange multipliers which should be put to zero after variation. The 
equations of motions for A^, cr, and pab (subject to detpa^ = 1) can be obtained by varying 
( p.23| ), with A" = 0. Here we present the results only; 

(a) D_(eVfei^+-) + e''{D_a){daa) - 9,(e'^D_a) - ^e>V(^-Pbd)(5aPce) 

+9fc(eV=Z}„pJ=0, (2.24a) 



[h] D+(e2>,,F^._^) - e'^{D+a){daa) + d^ie'^D^a) + le>V(^+P6d)(5aPee) 

-9fe(eV^D+p,e) =0, (2.24b) 

(c) {D+a){D_a) + 2D^+D_)a + ^p'''p'"'{D+pac){D^pM) + e" PabF^^F^l 

= 0, (2.24c) 

(d) D(+(e>"^D_)p,,) - ^e'^{pbcF^^F^^ - ^6%p,,F^^F^^) = 0, (2.24d) 

where the symmetric symbol is normalized such that {af3) := (a/3 + (3q)/2. Together 
with the 3 equations C± = 0, Cq = that we obtain by varying (|2.23|) with respect to 
A^, A°, these fleld equations are identical to the 10 Einstein's equations spelled out in the 
double null gauge, which we obtain by flrst varying the general E-H Lagrangian density 
( |A32|) in Appendix A, and then imposing the double null gauge ( p.ll| ). Therefore the 



Lagrangian density ( p.23[ ) is equivalent to the general E-H Lagrangian density ( |A32| ) , with 



the understanding that A°"s are to be set to zero after variation. 

The Lagrangian density ( p. 231 ) may be naturally interpreted as the Yang-Mills type 



Lagrangian density on the (l+l)-dimensional "space-time", interacting with the (1+1)- 
dimensional "matter" flelds a and pab- The corresponding local gauge symmetry is the 
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built-in diffA'2, and the "matter" fields couple to the diffA''2-valued gauge fields through 
the minimal couplings. In addition, each term in ( p.23| ), including the auxiliary equations, 
is manifestly invariant under the diffA^2 transformations. Therefore (|2.23| ) should be duly 
regarded as a gauge theory formulation of the vacuum general relativity 0-§l. 

That pab is the physical degrees of freedom can be also seen as follows. In this (1+1)- 
dimensional interpretation, the diffA^2-valued gauge fields A_^ are auxiliary fields since they 
have no propagating (i.e. no transverse traceless) degrees of freedom. Moreover, as we have 
seen already, a is also an auxiliary field that is determined by pab through the equations 
C± = 0. This confirms that the two physical degrees of freedom of gravitational field 
are indeed contained in pab- It seems appropriate to notice here that, in the propagating 
equations of motion (|2.24d| ) for pab, the source term is given by 

le^^PbcF^^F^^, (2.25) 

whose trace is precisely the local gravitational energy density, as we shall see in the next 
section. This indicates that the local energy density in general relativity indeed plays the 
analogous role as the local charge density does in Maxwell's theory. 

III. THE LOCAL GRAVITATIONAL ENERGY DENSITY 

In this section, we shall find the local Hamiltonian density of general relativity. This 
can be obtained simply by writing the local Lagrangian density (|2.23|) in the first-order 



form using a suitable time coordinate. The most natural time in this formulation seems the 
retarded time u [0,^. With the retarded time as our clock, the first term in (|2.23|) may 



be written as 

1 



-^YM :— -^e '^pabF^-F_ 



2a ^ jp a jp b 

1 



= e"" pabF^^[d+A^ - d-A^ - [A+^AJl'^ - -F^^). (3.1) 

In terms of the phase space variables (n^, A^), where Ua is defined as 
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Ha = e^^PahF^l, (3.2) 

this can be written as 

Cym = i^ad+A^ - V^v^^n.n, + A/D_n„ (3.3) 

ignoring the surface terms. Here -D-IIa is the diffA^2-covariant derivative of the density Tla 
defined as 

D_n, = 9_n,-[A_,n],, (3.4) 

where [A., n]a is the Lie derivative of Ha, 

[A_, n], = A^d^Iia + {daA^)Ii, + {d,A^)Iia. (3.5) 

The second and third term in ( ^.23| ) are aheady in the first-order form, apart from the 
terms proportional to A^ whose variation yields the Gauss-law equations associated with 
the residual diffA^2 invariance. Putting these all together, the Lagrangian density ( p.23|) can 
be written in the following Hamiltonian form^ 

C = U,d+A^ - e^(D_a)(9+cr) + ^e>V'(^-P6d)(5+Pac) 

-U-^-p-^IiJi, + A^Ca+ E ^"C'„, (3.6) 

^ a=±,0 

where Ca is given by 

Ca = D_Ua + e"{D_a){daa) - 9,(e-D_a) - ieVy^(D_p,rf)(9,Pee) 

+ab(eV^D_p,,), (3.7) 



which is the same as ( p.24a|) if we use (|3.2| ). From (|3.6| ) the local gravitational Hamiltonian 
density Ti is given by 



^Notice that the Hessian of this Lagrangian density is zero, so that the usual method of Haniil- 
tonization does not work. This is the reason that we did not introduce the momenta conjugate to 
a and pab- See for instance |36-^. 
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^ = ^e- VUn^ - A+^'a - E ^"Ca- (3-8) 

Using the 5 equations Ca = 0, C± = 0, and Cq = 0, the local Hamiltonian density ( p.8| ) 
becomes 

^ = ^e-^V^^n^n^ = ie2>„,F+_«F+_^ > 0, (3.9) 

which is positive-definite for any a and Ua, since the conformal 2-nietric Pab has a positive- 
definite signature. Thus, at least formally, we have obtained, in the double null gauge, a 
positive-definite local gravitational energy density for the vacuum general relativity! The 
time function associated with this non-zero local energy density is the retarded time u that 
we may choose at willQ. This is our proposal of the positive-definite local gravitational energy 
density in this article. This seems to be against the usual argument that, in general relativity, 
local gravitational energies can not be defined because they can be always "transformed" 
away due to the equivalence principle, let alone the positive-definiteness. In our definition 
of the local gravitational energy density, however, the "field strength" F,_^ in (^.13a|) is the 



coefficient of the commutator of the two null vector fields d± — A^da, which measures the 
twist of the parallelogram made of two successive parallel transports of these null vector fields 
along each other. Certainly the twist of this null parallelogram can not be "transformed" 
away even in a local Lorentz frame, and thus can serve as a measure of gravitational energy 
associated with the parallelogram surrounding the space-time point under consideration. 
The proposed local gravitational energy density is just the square of this local "field strength" 
multiplied by the canonical integration measure. 

In order to appreciate what this really means, however, we have to first define the volume 
integral of the local gravitational energy density S over a 3-dimensional hypersurface defined 



^For the asymptotically flat space-times, the number of possible choices of the retarded time u 
is equal to the number of an arbitrary, monotonically increasing function of three variables {u, y"') 



|11|. This may be true for other space-times as well. 
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by M = constant, and evaluate it for the asymptotically flat space-times, since the total 
gravitational energy is well-defined only for the asymptotically fiat space-times. As we now 
show, for the asymptotically fiat space-times, the volume integrals of the proposed local 
energy density over suitably chosen 3-dimensional hypersurfaces can be re-expressed as the 
Bondi and ADM surface integral at null and spatial infinity, respectively. We shall also 
derive the Bondi mass loss-formula as a negative-definite flux integral of a bilinear in the 
gravitational currents at null infinity. 

A. The Bondi Mass 

In this subsection we wish to show that, for the asymptotically fiat space-times, the 
volume integral of ( |3.9| ) over the u = constant null hypersurface is precisely the Bondi mass 
as measured at null infinity. Let us first notice that the volume integral E, where 

E=^J dvd'y e^'^pabF^^F^!^ > 0, (3.10) 



is positive-definite for any topology of N2 |T0[ . In order to express ( |3.10| ) as a surface integral, 
it is necessary to write it in a slightly different form using the field equations. For this let 
us consider the following identity 

D+D^a - D^D+(j = -F+^daCr - daF^^. (3.11) 

Using (|3.11|), the integral of the equation Cq = in (|2.2(Jb|) over N2 with the integration 



measure e" may be written as 
]^jd^y e'^'pabF^^F^' = J d^y e'^ji^a + 2{D^a){D_a) + 2D+D_a} + J d'y 9,(e'^F+_«). 

(3.12) 

The last term in (|3.12|) is zero for any 2-surface N2 that we assume compact without bound- 
ary. Thus (13.101 ) becomes 



E= f dvd^y e''[R2 + 2{D+a){D_a) +2D+D_a}. (3.13) 
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Let us also integrate the equation (|2.24c|) over N2, using (p.ll|) , to obtain 



j (fy e''[{D+(y){D_a) + 2D+D_a] 

-Jd'ydaie-F^^), (3.14) 



where the last term may be also dropped. Thus the volume integral ( p. 131) becomes 



E = Jdvd^y e'^{i?2 + {D+cj){D^a) - ]^p''''p'\D+pac){D.p, 

- J dvd^y e^ V,,F+_"F+_^ (3.15) 



b 

+- 



E = ^Jdvd^ye^''PabF^^F^ 

= ^-j dvd'^y e-^j/^s + {D^a){D_a) - ^pV'(^+Pac)(I^-Pfed)}. (3.16) 

To show that this[| can be expressed as a surface integral, the covariant null tetrad notation 
that we described in Appendix B is useful. Let us notice that the Gauss equation in the 
(2,2)-splitting of space-time is given by0 P0|-^ 

R2 + {D+a){D_a) - ]^p'^' p'\D^pa,){D_p,^) = h^'^^h^'^CABCD. (3.17) 

where Cabcd{.AiB,- ■ ■ = 0,1,2,3) is the conformal curvature tensor, and h^B is the 2- 
metric on the transverse surface A^2; i-e. the covariant form of (pat- Then the volume integral 
E may be written as 



^Notice that this energy integral is different from the one in Hayward's paper pH]. The energy 



density he proposed is the minus of Co in ( 2.19 ), modulo the Euler density, and is not positive- 



definite. This difference may be traced back to the fact that in his paper both u and v are treated 
as the time variables. 
^Here we used the vacuum Einstein's equations. 
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E = ^-jdv(fy e'^h^^h^'^'CABCD. (3.18) 

In Appendix C, we have shown that, using the Bianchi identity 

"^[mCabicd = 0, (3.19) 

this can be expressed as the surface integral 

E = ^l dv(fy e'^h^^h^'^'CABCD = ^ Hm vj (fy e''h^''h^''CABCD, (3.20) 

where hm means that the integral over N2 is to be evaluated at the limiting boundary 
value(s) of v. This expression picks up the coefficient of 1/f-term, and becomes precisely 
the Bondi mas^ |]40|-^ in the limit as v approaches to infinity! Notice that the parameter 
v/\/2 becomes the area radius in the limit f ^ cxd (keeping u = uq = constant), as our 
metric ( p^) approaches to 

ds^ -^ -2dudv + -{v- ufidd'^ + sin^M'/) (3.21) 

at null infinity. 

B. The Bondi Mass-Loss Formula 

In this subsection we shall continue to obtain the Bondi mass-loss formula in the presence 
of the gravitational radiation in the asymptotically fiat space-times. For this, we may 
simply take a w-derivative of the integral ( 3.2(J| ), using suitable vacuum Einstein's equations. 



However, since we wish to account for the mass-loss in terms of the physical degrees of 
freedom pab-, we shall work with the volume integral ( |3.13|) . 



Recalling that diffA^2 is the residual gauge symmetry of the metric ( p.2|) , we may fix this 
symmetry by choosing A^ = 0^ by a suitable coordinate transformation on N2. Then the 



^*^We have a factor of 1/3, which could be taken care of by a suitable normalization. 

^^This is equivalent to the assumption that d/du is a twist-free null vector field p|,p,44|. However, 



there could be some topological obstructions against globalizing this choice. 



equation (7+ = in ( |2.22|) reduces to the following Schrodinger equation 



S^e^/^ ^ ^2^.12 ^ Q^ ^j^g^g ^2 _ \p-^pC'^[Q^p^^)[Q^p^^) > 0. (3.22) 

o 

In the following we shall use this equation when we examine the rate of change in E as the 
retarded time u advances. Let us further notice that the metric (|2.2|) becomes, in the gauge 
A^ = 0, 

ds'^ = -2dudv + e^pabiA^dv + dy''){A^dv + dy^). (3.23) 

With y°' = {{}, ip), where '&, ip are the angles of 5*2, we find by comparing ( p.23| ) with ( p.21| ) 
the following asymptotic behaviors of the metric as v approaches to infinity, 

e'^ = 0(t;^), pab = 0{l), A^- = Oil/v^). (3.24) 

In the gauge A_^ = CQ, the volume integral ( p.l3| ) becomes 

E= f dvd^y e"{i?2 + 2{d+a){D_a) + 2d+D_a} 

= I dvd'^y[e''R2 + 2d+D_e''}, (3.25) 

where we used the identity 

a+D_e" = e"(9+a)(D_a) + e"9+D_a. (3.26) 

The first term in (|3.25| ) is the f-integration of the Euler number x, where 

X = ^Jd''y e''R2 = 2, 0, -2{g - 1), (3.27) 

for N2 = 5*2, T2, and the 2-surface S^ of genus g, respectively. Since the w-derivative of the 
Euler integral is zero, the rate of change in i^ as -u increases comes from the second term in 
(^.25|). For the asymptotically fiat space-times, we may assume N2 = S2 so that 



^^This gauge choice is only for convenience, since we already obtained the covariant expression of 
the Bondi mass in the previous subsection. 
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= 2 / Sy (aje") - 2 / £-y (S^e") (l + 0(«„-') 

Jv=oo,S2 Jv=vo,S2 



(3.28) 



,S2 ' Jv=V0,S2 

where the domain of the f- integration was chosen from wq to oo, and we used the asymptotic 



behaviors ( p. 241) . Here vq is some point that hes sufficiently far away from the sources of 
gravitational waves along the out-going null direction such that the gravitational waves are 
contained entirely in the range f o < f < oo at the instant u = constant. In this asymptotic 
region, we may also assume the out-going null condition [^,||,|Tl| , 






(3.29) 



for some function f[u, ^, if). From this, it is found that 



a^ e" = 2e"/252 g<x/2 (i + o(i/, 



--eW'(5+Pae)(9+p,,)(l + Oil/V) 



(3.30) 



where in the second line we used the Schrodinger equation (|3.22|) . Thus dE/du becomes 

-d^ = "2 iHS>i d"y eW^(5-fPae)(9+p,.) 



+ 



^ / d'y eW'(5+Pac)(5+p,,)(l + Oiv,')). 

/ Jv=vn.S9 



(3.31) 



/ Jv=Vo,S2 

Since there are no propagating gravitational degrees of freedom in the region v < vq the 
currents of gravitational waves must vanish in this region, so that 



ab 



p^'d+pac = for V <vo. 



(3.32) 



Thus the second term in ( p.31|) vanishes, and we finally have 

-^hm / AeVV'(5+Pac)(9+pb,)<0. 
2 ^-*°o Js2 



dE 
du 



(3.33) 



Apart from the integration measure e°", this fiux integral over 5*^ at null infinity is expressed 
entirely in terms of the physical degrees of freedom, and is negative-definite. This is precisely 
the Bondi mass-loss formula! It must be stressed that the gravitational energy carried 
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away to null infinity by the gravitational radiation is given in a bilinear combination of the 
gravitational currents, in excellent accordance with our experience that observables are very 
often expressed in bilinears of the physical fields. This strongly supports the view held by 
the geometric quantization school p6|-p^] that, in general relativity, as in other non-linear 
field theories, it is the current rather than the conformal 2-geometry that should be regarded 
as the fundamental physical field. 

The total radiated gravitational energy to null infinity between the null time interval mq 
and u can be obtained by integrating (|3.33|) , and is given by 

E{u) - E{uo) = -^limjyuj^^ (fy e>V'(5+Pac)(<9+pM). (3.34) 

C. The ADM Mass 

Now we shall show that the volume integral of the proposed local gravitational energy 
density ( p.lOp over a spacelike hypersurface reproduces the ADM surface integral. Let us 
make the following coordinate transformation 

1 
r = —-U + V. (3.35) 



In the new coordinates {u,r,y"'), the metric (2^) becomes 

ds'^ = -du'^ - 2dudr + e" pab[{A^ + -A^)du + A^dr + dy"] 

[{A^ + -A^)du + A^dr + dy^]. (3.36) 

Since we still have the residual gauge symmetry associated with the diffA^2 invariance in 



( p.36| ), we may well fix this residual symmetry by choosing 



A^ + ]-A^ = 0. (3.37) 

Then (|3.36|) reduces to 

ds"^ = -du^ - 2dudr + e'^p^^ (^A^dr + dy") {^A^dr + dy^) . (3.38) 
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In the limit as r -^ oo, both ( p.36| ) and (|3.38|) approach over to the flat space-time metric 



ds'^ -^ -du^ - 2dudr + -r\d^^ + sm^M<^^), (3.39) 

showing that r j \/2 becomes the area radius in this hmit. Moreover the w-coordinate is the 
proper time at each point on the u = constant hypersurface, suggesting that the volume 
integral be deflned over the spacelike hypersurface u = constant in the new coordinates, 
since the ADM surface integral is associated with a unit time translation at spatial inflnity. 
To flnd the relevant Hamiltonian for the ADM mass we need to write the local Lagrangian 
density in the new coordinates. The local Lagrangian density can be found directly from 
(S) 

^ ^ a=±,0 

(3.40) 

with 9_ replaced by dr everywhere. Thus the relevant volume integral is given by 

E = ^Jdrd^y e>,,F_,_«F+_^ |a_=a^ 

= ^j drd'^y e'^ji^a + {D+a){D^a) - lpV'(^+Pac)(/^-P6d)} \d.=d. 

= ^J drd% e''h^''h''''CABCD. (3.41) 

Repeating the same reasoning as in Appendix C, with Vt^^ = r, we flnd that the volume 
integral ( |3.41| ) becomes the surface integral, 

E = him rj d^y e'^h^'^h^'^CABCD, (3.42) 

which becomes precisely the covariant expression of the ADM mass of the asymptotically 
flat space-times in the limit as r approaches to inflnity! P0|-|i3 



IV. DISCUSSIONS 

In this article, we combined the double null formalism of general relativity with the KK 
formalism in the (2,2)-splitting, and proposed a local gravitational energy density of general 
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relativity. As we have seen so far, there are a number of notable features of this description 
which deserve further remarks. First of all, this formalism explicitly brings out the gauge 
theory aspects of general relativity of the 4-dimensional space-times. Although it has been 
realized for a long time that the local diffeomorphism invariance in general relativity is on 
a par with the local gauge symmetry in gauge theories, it seems fair to say that the full- 
fledged gauge theory formulation of general relativity is still lacking. Our 4-dimensional KK 
approach to general relativity in the (2,2)-splitting, using the double null gauge, seems to 
provide such a formulation, as we have described in this article. Thus we may well take care 
of the Gauss-law equations associated with the diffA^2 invariance by considering the diffA'2 
invariant quantities only. 

Moreover, this formalism shows that, in the double null gauge, local gravitational energy 
density of general relativity can be well-defined, and moreover, is positive-definite. The 
volume integral of this local gravitational energy density over the 3-dimensional null and 
spacelike hypersurface correctly reproduces the Bondi and ADM surface integral at null and 
spatial infinity, respectively. The Bondi mass-loss due to the gravitational radiation in the 
asymptotically fiat space-times is given by a negative-definite flux integral of the bilinear 
in the gravitational currents at null infinity. It should be mentioned that the proposed 
gravitational energy density can be also used to define quasi-local gravitational energies for 
a finite region of a given 3-dimensional hypersurface in a straightforward way [^ ,^5|-^ . 



The non-zero local Hamiltonian density proposed in this article also has a direct bearing 
to the problem of time []50|-p^ . The time associated with the non-zero Hamiltonian is the 



retarded time u, which has the physical meaning as the phase of the gravitational radiation 
when gravitational waves are present. The canonical analysis of our formalism is under 
progress |^, which will shed further light on this important issue. 

In addition, this formalism does seem to indicate the intriguing possibility that quantum 
general relativity of the 4-dimensional space-time may be regarded as a (l+l)-dimensional 
quantum field theory. For instance, one might even speculate that quantum gravity might 
be a finite theory, given that the renormalizability depends critically on the dimensions of 
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"space-time" . However, it must be addressed that this formahsm is for the vacuum general 
relativity only. It certainly is an interesting question to see whether this formalism can be 
extended to include matter fields. We leave this problem for the future investigation. 
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APPENDIX A: THE E-H LAGRANGIAN DENSITY IN THE GENERAL 

(2,2)-SPLITTING 

In this appendix, we shall make a general (2,2)-decomposition of space-time, and obtain 
the corresponding E-H Lagrangian density 0-^] without picking up a particular gauge. By 
examining the transformation properties of the metric in the (2,2)-splitting under the diffiV2 
transformations, we shall find that each field can be identified either as a scalar, a tensor, 
or a gauge field with respect to the diffA^2 transformations, respectively, suggesting that the 
KK program works even in the absence of any Killing vector fields. Then we simplify the 
general E-H Lagrangian density by introducing the double null gauge to obtain Cq in ( |2.12| ) 
0, which also has the diffA^2 symmetry as the residual symmetry. 

The 4-dimensional space-time may be regarded as a fibred manifold, i.e. a local product 
of two 2-dimensional submanifolds Mi+i x A^2, for which we introduce two pairs of the basis 
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vector fields (9^ = d/dx^{n = 0, 1) and da = d/dy°'{a = 2, 3), respectively. The corresponding 
metrics on Mi+i and A^2 will be denoted as •y^i, and (pab, respectively. Then the general line 
element of the 4-dimensional space-time can be written as 

ds^ = -f^^dx^'dx'' + (pabiA^^dxf" + dy''){Aj'dx'' + dy^). (Al) 

Formally this is quite similar to the "dimensional reduction" in KK theory, where Mi+i is 
regarded as the (l+l)-dimensional "space-time" and N2 as the "internal" fibre space. In 
the standard KK reduction one assumes a restriction on the metric, namely, an isometry 
condition, to make A^ a gauge field associated with the isometry group. Here, however, 
we do not assume such isometry conditions, and allow all the fields to depend on both x^ 
and y"". Nevertheless A^{x,y) can still be identified as a gauge field, but now associated 
with an infinite dimensional diffeomorphism group diffA'2. To show this, let us consider the 
following diffeomorphism of A''2, 

y'' = y\x,y), x^ = x^. (A2) 

Under this transformation, we find 

ip,u{x^y')=lpiu{x,y), (A3) 



QyC Qyd 
Qy'a Qy'. 



(t)'abi.x,y') = Tr77;Tr7i<Pcd{,x,y), (A4) 



A;,^{x,y') = ^-^A^^{x,y)-d,y\ (A5) 



such that the line element (|A1|) is invariant. Under the corresponding infinitesimal trans- 
formation 

5y'' = C{x,y), (5x^ = 0, (A6) 

these become 

81 ^.u = -[e, l^.u] = -Cdcl^^u, (AT) 

6(t>ab = -[e, <P]ab = -edc<Pab - {daC)<Pcb ' {d^e)<Pac, (A8) 

m; = -9^r + [A,, iT = -d,c + {A'^^d^c - r^c^;), (A9) 
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where the bracket represents the Lie derivative that acts on the "internal" indices a, b, etc, 
only. Notice that the Lie derivative, an infinite dimensional generalization of the finite 
dimensional matrix commutators, appears naturally. Associated with this diffA^2 transfor- 
mation, the diffA^2-covariant derivative D^^ is defined by 

D^ = d^-[A^, ], (AlO) 

where the bracket is again the Lie derivative along A^ = A^da- With this definition, we 
have 

SA; = ~D^C, (All) 

which clearly indicates that A"" is the gauge field valued in the infinite dimensional Lie 
algebra associated with diffA^2- Moreover the transformation properties (|A7|) and ( |A^ ) show 
that 7^jy and (pab are a scalar and a tensor field, respectively, under diffiV2. The field strength 
F^J" corresponding to A^ can now be defined as 

[D^, D,] = -F^fda = -{d,A: - d,A; - [A^, A,r}da, (A12) 

which transforms covariantly under the infinitesimal transformation ([AB|), 

5F,J' = -[i,F,^r. (A13) 

To obtain the E-H Lagrangian density, we have to first compute connections and curvature 
tensors. For this purpose it is convenient to introduce the following horizontal lift basis 



Qa = {df,, da) where |2 



d^ := d^ - A;da, da := da . (Al4) 

From the following commutation relations 

[dA,dB] = fAB''dc, (A15) 

we find the structure functions fAB^i^iV) 
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r a rp o, 

f ^ = —f ^ = ri A '' 

/^^^ = 0, otherwise. (A16) 



The virtue of this basis is that it brings the metric (|A1|) into a block diagonal form 



gAB=\ , (A17) 

V 0,6/ 

which drastically simplifies the computation of the scalar curvature. 

The connection coefficients and the curvature tensors in this basis are given by ||3^ , |56[1 

^AB = 2^ [pAgSD + dsgAD - dogAB) + ^g [IaBD - IbDA - IaDb) , 

p D _ f, -p D a r D j_r D-p E f^ D-p E . r E-p D 

^ABC — ^Bi^ AC ~ ^Ai^ BC + ^ BE ^ AC ~ ^ AE ^ BC + JAB ^ EC ■> 

Rac = g Rabcd, 

R = g^^RAc, (A18) 

where /abc '■= gcnfAB^- I^ components the connection coefficients are given by 

f ,." = \<t>'%<t>ac + \d,A; - ^0'"=0ae9eA;, 
f „; = l<P'%<Pac - \daAi - ^0^^0.e5.A;, 



\r''dA,, + \Y''Kcd,A: + ]^ 



l^^'dAa^ + -l^''<\>acd^.A: + -T'^'^bc^a^.^ 



Tafe' = ^0''('9a0M + 550,, - 9,0,,) , (A19) 



The following identities are also useful; 



1 ..;. ^ a 1 ... ^ „ 1 

1 



^^!^ = 7,l"^d.lap, Tpf = :.7"^5a7a/3, T a^ = l^^'dAa, ' d^A^ , 



^J = ^r%<Pab, rj = -r'd,<Pac. (A20) 
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The Ricci tensors and the scalar curvature are given by 

Rfiv = R^av + ^nau i ^ac = Rabc + ^aac i R = 'J R^u + (p Rac- (A21) 

Thus, in order to obtain the E-H Lagrangian density, we need to calculate 'j^'^R^u and (j)^'^Rac 
only. Let us first define R' and Rac as follows; 

p' f) V " f) V "-i-F ^F '^ F '^F f^ 

Rac = d,tj - dat,,' + f^.'^f,,'' - tjt,f. (A22) 

Notice that formally R' is identical to the Ricci tensor of Mi+i, except that d^ was used 
instead of d^. For this reason it might be called the "gauged" Ricci tensor of Mi+i, whereas 
Rac is the usual Ricci tensor of A^2- After a long computation we obtain '-f^^R^^v and (f)"''^Rac 
as follows; 

-Y''r\d,M{dcA:)+Y\daA;){d,A,') - ]p^^{daA;^){d,A^^) 
-\l'"'r'MdaA;){d,A^') 

-(V, + f ,;)(^7^>"'^.0,, - r^'daA:) 

-(Va + f,.")(^0'^V^5,7,.), (A23) 

r^'Rac = r^Rac + \r''r^<l>abF^c^F,^' - Ir'^'r^idal^a^dtlup) 
+ \r''Y''r^{dal,u){dblaf^) 

-(v^ + r^;){^Y''r%<Pat - Y'daAj^) 

-(Va + f,.")(^0V^9,7,.). (A24) 

Here the derivatives V^ and Va are compatible with the metrics 7^^ and 0bc in the horizontal 
lift basis, respectively, such that 

V^7a/3 = d^-fap - t'JlSp - t'JlaS = 0, 

^a(pbc = da(pbc - Tj(l)dc - Tj(j)M = 0. (A25) 
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For an object of the mixed-type such as X^.^.'^f'.'" in this basis, these derivatives act as follows; 






d^x^,':,^,,. 



-p 5 y C-/3-' 



+ TJX, 



Py C-5-- 



+ 



o y c---f3--- -p dy C---/3--- I -p cy d-/3 



ab d---a- 



fiS b---a 

ad b---Of- ' 



(A26) 



For instance, the followings are true. 



VM^'d^ct^ab) 



5^(r'40a6)-V(r'4^ 






(A27) 



which we used in (|A23 ) and ( A24 ). The scalar curvature R of the metric (^) becomes, 
using ([X23D and (|X2l , 






Here j^ and j" are defined as 



f = 7^>«''9,0,, - 27'^^^^^/, 
f = rVdb^,., 



and D,,(h„h is the diffXo-covariant derivative 



^fj.9ab 



Df,^ac = d^^ac - {daA^)^ec " {dcA^)^a 



d, 



[A 



}iy^ac [^ -^/i; V^Jacj 



where [A^, (p]ac is the Lie derivative of (pac along A^ = A^de 



[A 



fi-j V^Jac 



^/5e0ac+(5a^/)0ec+(5A 



(A28) 



(A29) 



(A30) 



(A31) 
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Thus the E-H Lagrangian density in this (2,2)-sphtting finally becomes 



C' = ^^^fpR 



-^J(j) 



r 



flU 



RL + r'Rac - T7^'^7"^<^a6i^„ "i^, 



fihi 



fia v(5 






(A32) 



where 7 = det7^i, and = det^^fe. It can be shown that the last two terms in (|A32|) are 
total divergences P-§|, 

(A33) 
(A34) 



-7v/0(v, + f,;)/ = a,(v/^v^/) - da[V^4^A;f), 



using ( [A20| ) and the following identities 



v,f = d,f + r^^f 



afj. J 5 
a I -p b -a 



Var = dar + Tjf. 



In the {u, v) coordinates 



1 



U = —;=(x^ — X^), 

V2 



x/2 



'x'^ + x'] 



the following substitution 



7+- = -1, 7++ = 7- 



0, 



together with 



a^ = -^{a,^ta,'^), 



(A35) 



(A36) 



(A37) 



(A38) 



leads to the double null gauge ( |2.2| ) and enormously simplifies the E-H Lagrangian density 
( |A32| ). The resulting expression is Cq in ( p.l2| ). 



30 



APPENDIX B: THE COVARIANT NULL TETRADS 

In this appendix, we describe the kinematics of a Lorentzian space-time of 4-dimensions 
using the covariant null tetrads 0-|T^|T^. This allows us to compare the variables in the 
traditional double null hypersurface formulation and our KK variables directly. Moreover, 
in order to express the volume integral ( ^.181 ) as a surface integral, it is better to use the 



covariant null tetrad notation. Let the two real dual null tetrads I a and ua (A=0, 1,2,3) be 
the gradient fields for some scalar functions u and v, 

Ia = ^au, Ua = Vav, (B1) 

so that V[bIa] = ^[BnA] = 0. The dual vector fields du and dv are related to the dual null 
tetrads by 

du = lAdX^, dv = UAdX^. (B2) 

We also have the vector fields d/du, d/dv, and d/dy"' {a = 2, 3) which we may write 

(^ _ A ^ 9 _ A ^ 9 _ A 9 (Y>n\ 

If we choose the basis vector fields of space-time so that Va = {d/du,d/dv,d/dy°'), the 
components of the dual null tetrads and vector fields are given by 

/a = (1,0,0,0) nA = (0,1,0,0), 

«^ = (1,0,0,0), t;^ = (0,1,0,0), l// = 5/. (B4) 

From this it follows that 

l^u^ = haV^ = 1, IaV^ = uaW^ = 0. (B5) 

Since IaI^ = uau"^ = 0, the null tetrad l^ and n^ may be chosen as 

ax^ dv - dy"' dX^ du + dy"' ^ ^ 



I.e. 
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/^ = (0,-l,A_«), n^ = (-l,0,A/), (B7) 

such that 

/^n^ = -1. (B8) 

[n and / are the minus of the horizontal hft vector fields d± = d± — A^da in the [u, v)- 
coordinates, respectively.) The condition (|B^ ) is equivalent to the previous normalization 
condition g^ = —1 in ( p.3|) , and means that, given an arbitrary function u, the function 
V must be chosen in such a way that the normalization condition ( pSD is satisfied. We still 
have the freedom to orient l^ and n-^ in space-time, and we fix this freedom by demanding 
that l^ and ra^ are normal to the 2-dimensional spacelike surface N2 whose tangent vector 
fields are d/dy"", 

hABl"" = hABu'' = 0, (B9) 

where Hab is the metric on N2 {Jiab is the covariant expression of (pab)- Using Kab and Ia-, 
Ha, the space-time metric qab may be written as 

QAB = flAB - {lATiB + UaIb)- (BIO) 

From these relations, we easily find that (|B10| ) is identical to the metric (pl^ ) 

d.s^ = -2dudv + (l)abiA_^du + A^dv + dy"-){A^du + A^dv + dy^). (Bll) 

APPENDIX C: THE BONDI SURFACE INTEGRAL 



We now show that the volume integral ( |3.18|) 



E = \j dvd^y e''h^''h''''CABCD (CI) 

can be expressed as a surface integral, using the Bianchi identity of the conformal curvature 
tensor 
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V[AfC^B]CD = 0. (C2) 

Let us notice that, due to the Bianchi identity, the following is true for any scalar function 

V[m{P~ Cab]cd) = (V[jv/fi~ jCab]cd- (C3) 

If we contract (C3) by h'^'-^h^'^, it becomes 



Let us choose Q~^ as a function of {u,v) only, and let Vm = V-. Then the r.h.s. of (|C4|) 
becomes 

since h^^V b^'^ = 0. The l.h.s. of (|C4D becomes 

h^''h^''V[-{n~'CAB]CD) = Vl-^h^'^h^^n-'CAB^CD) 

-V[^{h^''h^'')Q-'CAB]CD. (C6) 

Using h^'^' = g^^ + l^nP + n^l^ and the properties of the conformal curvature tensor 

9 ^ABCD — 9 ^ABCD — U, 

CaBCD = C[AB][CD], (C7) 



the second term in the r.h.s. of ( |C6D may be written as 

Vi^{h^^h'''')n-'CAB]CD = -2Vi^{lAnB]lcnD)n-'C^''^''. (C8) 

Since I a and ha are non-zero only when A, B are + or — , we have 

V[_(/AnB]/cnD) =0, (C9) 

due to the repeated indices in the anti-symmetric symbol. Therefore the l.h.s. of (|C4D 



becomes 
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/l^^/l^^ V[_ {n-'CAB]CD) = V[- {h^''h^''n-'CAB]CD) . (CIO) 

Thus (|C|) becomes, using (^) and (|CTO| ), 



+VB(h''''h^''n-^C.ACD). (Cll) 

Integrating ( |C11| ) over the u = constant hypersurface with the canonical measure \/h ||35|| , 
it becomes, using y/g = \fh for the metric ( P10|) , 



f dvd^yVhh^^h^^(V-n-^)CABCD = I dvd^yVhV-ih^^h^^n-^CABco) 

= hm / d^yy/h (^h^^h^^Q-^CABCD) , (C12) 



since the surface integrals coming from the last two terms in the r.h.s. of (|C11|) vanish for 



any 2-surface N2 compact without boundary. Here lim means that the integral over N2 must 
be evaluated at the limiting boundary value(s) of v. Now let us choose Q such that Q~^ = v. 
Then the identity ( pi2| ) becomes 



/ dvd^yVhh^^h^^CABCD = lim vf d^yVhh^^h^^CABCD- (C13) 

Since yh = e°" in our previous notation, the volume integral ( |3.18| ) becomes 



E=^J dvd% e''h^''h''''CABCD = I lim vj d^y e''h^^h''''CABCD, (CM) 

which is the desired expression for the Bondi surface integral as v approaches to infinity. 

APPENDIX D: THE POSITIVE-DEFINITENESS OF k^ 

Here we show that k^ is positive-definite. Let us introduce super indices A', B' for the 
symmetric combinations [ac] and ihd), respectively, so that 

pA' ■= Pac, PB' ■= Pbd- (Dl) 

Define the supermetric Ga'b' and its inverse G^ ^ by 
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Ga'B' := \{Pa,Pcd + PadPcb). G^' ^' .= ^(pV' + P V), (D2) 

such that 

G^'^^'Ge'b' = 5^B', where 5^, = ^{6\6^, + 5\5\). (D3) 

This supermetric raises and lowers the super indices 

G^''''pB'=p^\ Ga'B'p''' =PA'^ (D4) 

and has a positive-definite signature since it becomes 

G^'^' = diag(+l, +1/2, +1) for p"^ = (5"^ (D5) 

Therefore it follows that 

4 = lpV'(^±Pac)P±Pfed) 
o 

= Ig^''''{D^pa'){D^Pb') > 0. (D6) 
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